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Geometric frustration arises when lattice structure pre- 
vents simultaneous minimization of local interactions. 
It leads to highly degenerate ground states and, sub- 
sequently, complex phases of matter such as water ice, 
spin ice and frustrated magnetic materials. Here we re- 
port a simple geometrically frustrated system composed 
of closely packed colloidal spheres confined between par- 
allel walls. Diameter-tunable microgel spheres are self- 
assembled into a buckled triangular lattice with either up 
or down displacements analogous to an antiferromagnetic 
Ising model on a triangular lattice. Experiment and the- 
ory reveal single-particle dynamics governed by in-plane 
lattice distortions that partially relieve frustration and 
produce ground-states with zigzagging stripes and subex- 
tensive entropy, rather than the more random configura- 
tions and extensive entropy of the antiferromagnetic Ising 
model. This tunable soft matter system provides an un- 
charted arena in which the dynamics of frustration, ther- 
mal excitations and defects can be directly visualized. 

Geometric frustration arises in physical and biological 
systems' ranging from water^ and spin ice-' to magnets'*'^, 
ceramics*", and high-Tc superconductors^. The essence of this 
phenomenon is best captured in the model of Ising spins ar- 
ranged on a two-dimensional (2D) triangular lattice and in- 
teracting anti-ferromagnetically **•''; two of the three spins on 
any triangular plaquette within this lattice can be antiparallel 
to minimize their anti-ferromagnetic (AF) interaction energy, 
but the third spin is frustrated because it cannot be simultane- 
ously antiparallel to both neighbouring spins (Fig.[T]\). Such 
frustration leads to materials with many degenerate ground 
states and extensive entropy proportional to the number of 
particles in the system. Consequently, small perturbations can 
introduce giant fluctuations with peculiar dynamics. Tradi- 
tionally, these phenomena have been explored in atomic ma- 
terials by ensemble averaging techniques such as neutron and 
X-ray scattering, muon spin rotation, nuclear magnetic reso- 
nance, and heat capacity and susceptibility measurements^. 
More recently, artificial arrays of mesoscopic constituents 
have been fabricated in order to probe geometric frustra- 
tion at the single- 'particle' level. Examples of the latter in- 
clude Josephson junctions'", superconducting rings", ferro- 
magnetic islands'^''-''''*, and recent simulations'^ of charged 
colloids in optical traps. Observations in these model sys- 
tems, however, have been limited to the static patterns into 
which these systems freeze when cooled. Thus many ques- 
tions about frustrated systems remain unexplored, particularly 
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FIG. 1 : Ising ground state. (A) Three spins on a triangular plaque- 
tte cannot simultaneously satisfy all AF interactions. (B) For colloids 
confined between walls separated of order 1 .5 sphere diameters (side 
view), particles move to opposite walls in order to maximize free vol- 
ume. (C,D) Ising ground state configurations wherein each triangular 
plaquette has two satisfied bonds and one frustrated bond. (C) Zigzag 
stripes generated by stacking rows of alternating up/down particles 
with random sidewise shifts; all particles have exactly 2 frustrated 
neighbours. (D) Particles in disordered configurations have 0, 1,2, 
or 3 frustrated neighbours (red hexagons). 



those associated with single-particle dynamics. For exam- 
ple, how, when, and why do individual particles change states 
to accommodate their local environments, and what kinetic 
mechanisms govern transitions to glassy phases? 

Here we report on the static and dynamic properties of a 
self-assembled colloidal system analogous to Wannier's AF 
Ising model**. Densely packed spheres between parallel walls 
form an in-plane triangular lattice with out-of-plane up and 
down bucklin g'^-'^i'^i'^'^"'^''^^-^^'^"''^^'^^ . The up-down states 
of the spheres produced by buckUng are analogous to up-down 
states of Ising spins (Fig. [TJi). Nearest-neighbour excluded 
volume interactions between particles favour opposite states 
for neighbouring particles, as do the AF interactions between 
neighbouring spins in the Ising model. In contrast to engi- 
neered mesoscopic system s "''"•'^''^''"' , however, the colloidal 
system facilitates easy tuning of the effective AF interaction 
through changes in the diameter of temperature-sensitive mi- 
crogel spheres^^. The colloidal system also permits direct vi- 
sualization of thermal motion at the single-particle level. In 
the limit of weak confinement, or weak interaction strength, 
system properties closely follow those predicted for the AF 
Ising model, but in the limit of strong confinement, they do 
not. For strong interactions, the lattice deforms to maximize 
free volume, and the collective nature of the free-volume- 
dominated free energy characteristic of most soft-matter sys- 
tems becomes important. We understand these effects the- 




FIG. 2: Buckled 
monolayer of col- 
loidal spheres. 
Movies are in Sup- 
plementary Informa- 
tion. (32 /im)^ area 
at T = 24.7°C (A-C) 
and 27.1°C (D-F). (A, 
D): Bright spheres: up; 
: dark spheres: down. 

;; (B, E): Labyrinth 
patterns obtained by 
drawing only the frus- 
trated up-up (red) and 
down-down (green) 

,J bonds. (C, F): Cor- 

Sx responding Delaunay 
triangulations. Blue 
dots mark defects in 
the triangular lattice, 
i.e. particles that do 
not have exactly six 
nearest neighbours. 

"^I^ Thermally excited 
triangles with three 
'- spheres up/down are 
labelled by red/green. 



oretically in terms of tiling of the plane by isosceles trian- 
gles. The tiling scheme identifies a ground-state consisting of 
zigzagging stripes with sub-extensive entropy. Interestingly, 
in contrast to Ising-model predictions, first measurements of 
single-particle 'spin-flipping' suggest that flipping dynamics 
depend not only on the number of nearest-neighbour frus- 
trated 'bonds', but on how these bonds are arranged. Thus the 
paper begins to explore connections between frustrated soft 
matter and hard materials such as frustrated AF media. 

Experimental System. For walls separated by distances 
on order 1 .5 sphere diameters, the particles maintain in-plane 
triangular order (see Supplementary Information) but buckle 
out-of-plane (Fig. |2l A, D). This buckling minimizes sys- 
tem free energy, F — U — TS, where U is the internal 
energy, T temperature, and S entropy; spheres move apart 
to lower their repulsive interaction potential energy U and 
to increase their free volume V, which in turn leads to an 
entropy increase with S oc liiV. The effective repulsion 
causes spheres to move to the top or bottom wall, and nearest- 
neighbours maximize free volume by moving to opposite 
walls (Fig. [TJ}). Buckled colloidal monolayers were first ob- 
served more than two decades agc '^'^^i'^ , and the AF analogy 
was then suggested'^-^*^. However, to date, few quantitative 
measurements have been performed on this system class, and 
the themes explored by most of the early work centred largely 
on structural transitions exhibited by colloidal thin films as a 
function of increasing sample thickness ^^'^^-'^-^"'^^'^-^'^"^ , rather 
than their connection to frustrated anti-ferromagnets. The use 



of temperature-sensitive diameter-tunable NIPA (N-isopropyl 
acrylamide) microgel spheres^^ also distinguishes our experi- 
ments from earlier work. By varying temperature we change 
particle size and sample volume fraction and, therefore, vary 
the strength of the effective AF interparticle interactions. 

Samples were annealed at low volume fraction near the 
melting point to produce 2D crystal domains with ^lO'^ 
spheres covering an area of order (60/^r7i)^. Video microscopy 
measurements were carried out far from grain boundaries on 
a ^^(32/1771)^ central area (^2600 spheres) within the larger 
crystal domain. Particle motions were observed by micro- 
scope, recorded to videotape using a CCD camera and tracked 
by standard image-processing techniques'^. Our colloids are 
very weakly charged'^ and were measured to have short- 
ranged repulsive interactions^*^. Furthermore, NIPA spheres 
are nearly density matched in water, so gravitational effects 
are negligible. In most colloid experiments the important ther- 
modynamic control variable is particle volume fraction. The 
present experiment achieved substantial variation in sphere 
diameter using small changes in temperature, which altered 
thermal energies by less than 1%. In this paper we monitor 
and report temperature rather than volume fraction because 
the interactions between spheres contain a soft tail that intro- 
duces some ambiguity into the assignment of a geometric di- 
ameter to the particles. Below 24°C the system is jammed and 
no dynamics are observed. Above 27.5°C the in-plane crys- 
tals melt. Our primary measurements of the frustrated states 
probe five temperatures from 24.7°C to 27. 1 °C in 0.6°C steps. 



In this range, the hydrodynamic diameter of the particles de- 
creases linearly with increasing temperature from 0.89 /im to 
0.76 /im, while the average in-plane particle separation re- 
mains 0.7 fim. To reach thermal equilibration, the sample was 
annealed near the melting point before the temperature was 
slowly decreased. No hysteresis was observed when slowly 
cycling through this temperature range. 

Anti-Ferromagnetic Order. The images in Fig. |2] A, D 
show roughly half of the spheres as bright because they are in 
the focal plane of the microscope; the other half, located near 
the bottom plate, are slightly out-of-focus and appear dark. 
We discretize the continuous brightness profile of the particles 
into two Tsing' states with Si = ±1. The nature of the frus- 
trated states can be exhibited in different ways in processed 
images. One way focuses on the 'bonds' between particles. 
We refer to pairs of neighbouring particles in opposite states 
{SiSj = —1) as satisfied bonds, i.e. satisfying the effective AF 
interaction, and to up-up or down-down pairs (with SiSj = 1) 
as frustrated bonds. Images of these bonds show that the frus- 
trated bonds form a nearly single-line labyrinth (Fig. |3J) at 
low temperature that then nucleates into domains (Fig. |2ji) 
at high temperature. Local AF order is alternatively charac- 
terized by the average number of frustrated bonds per parti- 
cle, (Nf). In the limit of weak interactions, an Ising system 
chooses a completely random configuration with half of the 
six bonds satisfied and half frustrated, leading to (Nf) = 3. In 
the limit of strong interactions, on the other hand, each trian- 
gular plaquette has one frustrated bond (Fig.[TJ\), a third of the 
bonds are frustrated, and (Nf) = 2. (Nf) is a linear rescaling 
of the density of excited triangles (3 up or 3 down) in Fig. |2l 
C, F, which ranges from in the Ising ground state to 0.5 for 
a random configuration. We find that (Nf) decreased from 
approximately 2.5 to 2.1 in the temperature interval 27.1°C- 
24.7°C. Detailed statistics of the different local configurations 
are presented in Supplementary Table S 1 . 

We first consider the static properties of the frustrated sam- 
ples. In particular we aim to identify similarities and differ- 
ences between the colloidal system and the Ising model. As 
the temperature is lowered to increase the particle diameter, 
(Nf) is observed to approach 2. This behaviour is expected 
in the Ising model ground state. However, the vast majority 
of Ising ground-state configurations are disordered. The col- 
loidal monolayers, by contrast, condense into stripe phases. 
The stripes are not straight, as could be produced by higher- 
order interparticle interactions'". Rather they bend and form 
zigzag patterns?2i2i2i2^i2^ (see Fig.|2]\ and Supplementary Ta- 
ble SI). In this colloidal zigzag striped phase, we measured 
spatial correlations r{i — j) = [{siSj) — (s)'^]/[(s^) — (s)'^] 
over separations \i — j|, along the principal lattice directions, 
of up to 20 particles and found that they decay exponentially 
in magnitude with alternating sign (Supplementary Fig.S4). 
T{i — j) is positive for i — j even and negative for i — j odd. 
In contrast, T{i — j) averaged over the Ising ground state is 
positive when i — j is an integer multiple of 3. Furthermore, 
for zigzagging stripes each particle has exactly two frustrated 
neighbours (Fig. [if!), whereas in the fully disordered Ising 
ground-state Nf can be 0, 1, 2, or 3 (Fig.[TJ)) and only the 
average (Nf) is 2. These observations suggest that fluctua- 
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FIG. 3: Fluctuation in the number of frustrated bonds per par- 
ticle as a function of its average. Experiments quantitatively agree 
with hard-sphere simulations at different plate separations h, normal- 
ized by the average in-plane lattice constant L. Simulations collapse 
onto a single curve and deviate significantly from the behaviour in 
the Ising model. 



tions in Nf, i.e. Var(iV/) = (iV|) - (A^/)^, might be a useful 
measure for distinguishing the zigzag stripe phase observed 
here from the disordered Ising ground-state. Figure [3] plots 
the behaviour of Var(iVy) as a function of (Nf) for the Ising 
model and for data obtained both from experiments and from 
hard-sphere Monte Carlo (MC) simulations. Results from ex- 
periment and simulation agree at both low and high volume 
fraction and differ from those of the Ising model, especially at 
high volume fraction wherein interactions are strong. 

Zigzagging Stripes. Ideal geometrically frustrated sys- 
tems, such as the AF Ising model, are highly degenerate with 
extensive entropy at zero temperature. However, in real ma- 
terials, subtle effects, for example anisotropic interactions^ 
long-range interactions^', boundary conditions'^ and lattice 
distortions''^''"' relieve frustration. Our partially ordered 
zigzag stripe phase at high volume fraction is an example of 
frustration relief by lattice distortion. In the colloidal mono- 
layer the triangular packing is self-assembled, and (like atoms 
in real solids) the particles are not forced to remain at fixed po- 
sitions on the lattice^^. This deformability and the fact that the 
free volume of the system is a collective function of all parti- 
cle positions breaks the mapping to simple Ising models with 
pair-wise additive nearest neighbour interactions. In fact, the 
positions of the colloidal particles may be thought of as com- 
prising a planar structure that crumples between the two con- 
fining planes. This "crumpling" leads to deformations of the 
planar triangular lattice with satisfied bonds (projected onto 
the plane) on average 3 — 4% shorter than frustrated bonds. 
This difference is consistent with the notion that each pair 
of neighbouring particles prefers to be separated by the same 
fixed distance in 3D, whether or not their connecting bond is 
satisfied. 

A simple tiling argument demonstrates why the colloidal 
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FIG. 4: Tiling the plane with isosceles triangles. (A) Close-packed 
spheres are separated by one particle diameter d in 3D. This dis- 
tance projected on the 2D plane remains d for a frustrated bond 
(siSj = 1), but is reduced to a:: = y^d'^ — (/i — d)^ for a satis- 
fied bond (SiSj — —1). (B) Viewed from above, each plaquette in 
the lattice tends to deform to an isosceles triangle with one long side 
(d) along the frustrated bond and two short sides (x < d) along the 
satisfied bonds. The angle larger than vr/S is marked in red. (C) All 
possible in-plane local particle configurations appearing in the Ising 
ground state. The isosceles triangles can tile the plane without extra 
space only for Nf = 2. The "white space" for Nf = 0,1,3 corre- 
sponds to additional excluded volume. (D, E) Tilings corresponding 
to striped and disordered Ising ground-state configurations, respec- 
tively, of Fig. [T] C, D. 



system ground-state configurations of stripes and zigzags 
pack better than the disordered Ising configurations. Further- 
more, the tiling model shows explicitly that maximal volume 
fractions of stripe and zigzag phases are the same. Each tri- 
angular plaquette in the Ising ground-state contains two satis- 
fied bonds and one frustrated bond. Thus, when spheres are 
close-packed in 3D, the equilateral triangle defined by each 
such triplet of neighbouring particles is tilted, and when pro- 
jected onto the 2D plane, it deforms into an isosceles triangle 
with two short sides along the satisfied bonds and one long 
side along the frustrated bond (Fig. |4] A, B). Subsequently, 
close-packed configurations of the buckled spheres in 3D are 
described by tilings of the plane by isosceles triangles. Fig- 
ure|4p shows the configurations of isosceles triangles for dif- 
ferent numbers of frustrated bonds (iV/) in the basic hexago- 
nal cell. By summing up the angles around the central vertex, 
one immediately sees that for iV/ = 0,1,3, the triangles can- 
not close-pack. Only the two configurations with iV/ — 2 en- 
able tiling the plane with isosceles triangles, or, equivalently, 
close-packing of the buckled spheres in 3D. Configuration 2b 
corresponds to a bend in a stripe, and 2c to a stripe continu- 
ing along a straight line. Both have the same maximal volume 
fraction, thus corroborating observations of zigzagging stripes 
in the experiments and simulations. 

Experiments and simulations indicate a preference of the 
stripes to form straight segments rather than to bend eas- 



ily and thus to generate randomly zigzagging configurations 
(Fig.|2j\). Zigzagging stripes can be viewed as a random stack 
of ordered lines of alternating up and down particles (Fig. [if!), 
thus straight and zigzagging stripes are analogous to the face- 
centred cubic (FCC) lattice and the random hexagonal-close- 
packed (RHCP) structure's in 3D. Straight and zigzagging 
stripes are equivalent in the close-packed limit by having the 
same maximal volume fraction. However, for smaller volume 
fractions there may be an order-by-disorder effect^'^^, giving 
a small free volume advantage of straight stripes over zigzag- 
ging ones, similar to the free volume advantage^^ of FCC over 
RHCP in 3D. 

Instead of an extensive entropy at zero temperature-, 
wherein S scales linearly with the number TV of particles in 
the system, here the entropy is subextensive. The number 
of zigzagging striped configurations grows exponentially with 
the linear dimension of the system (there are two possible 
ways of placing one row relative to its predecessor in Fig.[T]Z!), 



hence the entropy scales as viV. Alternatively put, a non- 
branching single-line labyrinth is dictated by the particles on 
the boundary, and for the system to rearrange from one zigzag 
stripe configuration to another, a percolating cluster of order 
ypN particles should be flipped. 

Dynamics. Taken together these observations have inter- 
esting consequences for the ground-state dynamics of frus- 
trated systems. The Ising ground-state has a local zero-energy 
mode, as shown in configuration 3c in Fig.|5]\: the central par- 
ticle can flip without changing the energy of the system, thus 
rapidly relaxing spin correlations via a sequence of such sin- 
gle spin flips, even at zero temperature. For buckled spheres, 
on the other hand, the close-packed configurations have only 
particles with A^/ — 2, and, moreover, even a particle with 
N f = 3 in an excited configuration has to cross an energy bar- 
rier in order to flip. Like the glassy behavior of an Ising model 
on a deformable lattice-'^-''^, the slow dynamics we observe at 
low temperature is a consequence of the absence of local zero- 
energy modes in the bulk. Subextensive ground-state entropy 
also appears in related models emulating systems with glassy 
dynamics'*". 

Online movies (see Supplementary Information) permit us 
to directly visualize 'spin flipping' as well as the motions 
of thermal excitations and defects in frustrated systems for 
the first time. Thermal excitations labelled as coloured tri- 
angles in Fig. |2] C, F were typically found to be gener- 
ated/annihilated in pairs due to the flipping of a particle shared 
by the two triangles. Well-isolated thermal excitations, on the 
other hand, appear to be more stable. To quantify these effects, 
we first extract the full time trajectory, Si(i), of each particle 
i from the movies. In Fig. |5|J we plot the single particle au- 
tocorrelation function C(t) = [(sj(t)si(O)) - (si)^]/[(sf) - 
(si)^], averaged over all particles not at lattice defects. As 
the temperature is lowered, the correlation function develops 
a stretched exponential form, C{t) — exp[— (t/r)'']. The 
measured relaxation time t exhibits a dramatic increase as the 
particles swell at low temperature, while the extracted stretch- 
ing exponent /3 decreases, indicating slow dynamics similar to 
those found in glasses. 

To further explore the dynamics of different local configu- 
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FIG. 5: Single particle dynamics. (A) Local configurations are 
labelled by their value of Nf and an index a,b,c indicating the 
precise geometrical arrangement of the frustrated neighbours for 
Nf = 2, 3, 4. Symmetry to rotation and inversion reduces the 2 
possible configurations to the 13 given here. (B) Single particle au- 
tocorrelation functions plotted versus decay time. Lines are fits to 
stretched exponentials C{t) — exp[— (f/r)''], with r and (3 given 
in the inset. (C) Flipping rates for the different local environments. 
Configurations 7n and 5n are defects in the in-plane lattice, with 7 
and 5 nearest neighbours. 



rations (defined in Fig. EK), Fig. [Sf shows the flipping rate 
fr of single particles with a fixed neighbour structure. We 
measured the probability p that a particle flips between con- 
secutive images given that the Ising states of its neighbours 
remained unchanged. The time intervals of dt — 1/30 sec 
between frames were short enough such that p was typically 



small (0.36 at most) and the flip rate could be approximated 
by /,. — p/dt. At high temperature, the behaviour is simi- 
lar to that of an Ising model undergoing Glauber dynamics: 
/,■ ^ g-AB/fceT yyjjgj-g tjjg energy difference AE is propor- 
tional to the difference in Nf before and after flipping. In 
this regime, particles with large Nf flip slightly more slowly 
at higher temperature because of the weaker interactions be- 
tween spheres. As the volume fraction is increased by lower- 
ing the temperature, the particle dynamics slow by 1-2 orders 
of magnitude and, more interestingly, significant differences 
develop between different geometrical configurations with the 
same Nf .Such phenomena may not appear in the simple Ising 
model where the Hamiltonian depends only on Nf. 

Defects in the underlying lattice can strongly affect the 
properties of frustrated systems. However, detailed knowl- 
edge about the role of defects in frustrated systems is very lim- 
ited. Our experiments permit us to directly visualize defects 
nucleating, annihilating, and diffusing (see Supplementary In- 
formation movies) . By comparing trajectories containing dif- 
ferent numbers and types of defects, our initial studies suggest 
that defect particles have enhanced in-plane diffusion (Supple- 
mentary Fig.SS) and slower flipping dynamics than averaging 
over particles with six nearest neighbours. 

Conclusion. We have presented experimental measure- 
ments of single-particle dynamics in a geometrically frus- 
trated system. Colloidal spheres with tunable diameter self- 
assemble to buckled monolayer crystals and form a system 
analogous to the triangular lattice AF Ising model. By tuning 
the volume fraction, we found that at high compaction, in- 
plane lattice deformation relieves most frustration and yields 
a zigzag stripe ground-state with subextensive entropy. The 
'free spins' in the Ising ground state are removed; thus the 
system become glassy as the volume fraction is increased. A 
theoretical analysis shows that these features can be captured 
by a hard-sphere model. We measured spatial correlations and 
the statistics of various local configurations as well as their 
flipping rates and found strong dependences on arrangements 
of neighbouring particles. As the glassy phase is approached, 
we observed dramatic slowing of the dynamics and formation 
of stretched exponential correlation functions. Single-defect 
dynamics were directly visualized and measured for the first 
time. Defects have faster in-plane diffusion and slower out- 
of-plane flipping than the average. 

Our demonstration and analysis of this self-organized col- 
loidal 'antiferromagnet' opens the door for the study of de- 
tailed single-particle dynamics in frustrated systems and be- 
gins an exploration of the connections between frustrated soft 
materials and the more studied frustrated magnetic materi- 
als. Many further manipulations can be readily applied to 
this system, e.g. external control of particle motions with 
optical tweezers, gravitational fields, electric fields, different 
particle interactions, and defect doping. Theoretically, it will 
be interesting to consider possible modifications to the Ising 
model that generate a zigzagged stripe ground state and to 
study the glassy dynamics arising from its subextensive zero- 
temperature entropy. 
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